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Abstract 


Aircraft and ground-based radar data from the Tropical Ocean and Global Atmosphere Coupled- 
Ocean Atmosphere Response Experiment (TOGA COARE) show that convective systems are 
not always vertical. Instead, many are tilted from vertical. Satellite passive microwave 
radiometers observe the atmosphere at a viewing angle. For example, the Special Sensor 
Microwave/Imager (SSM/I) on Defense Meteorological Satellite Program (DMSP) satellites and 
the Tropical Rainfall Measurement Mission (TRMM) Microwave Imager (TMI) on the TRMM 
satellite have an incident angle of about 50°. Thus, the brightness temperature measured from one 
direction of tilt may be different than that viewed from the opposite direction due to the different 
optical depth. This paper presents the investigation of passive microwave brightness 
temperatures of tilted convective systems. 

To account for the effect of tilt, a 3-D backward Monte Carlo radiative transfer model has been 
applied to a simple tilted cloud model and a dynamically evolving cloud model to derive the 
brightness temperature. The radiative transfer results indicate that brightness temperature varies 
when the viewing angle changes because of the different optical depth. The tilt increases the 
displacements between high 19 GHz brightness temperature (Tbi 9 ) due to liquid emission from 
lower level of cloud and the low 85 GHz brightness temperature (Tbss) due to ice scattering from 
upper level of cloud. As the resolution degrades, the difference of brightness temperature due to 
the change of viewing angle decreases dramatically. The dislocation between Tbi 9 and Tbss, 
however, remains prominent. 

The successful launch and operation of the TRMM satellite provide us an opportunity to 
examine tilted convective systems using collocated radar and radiometer data. TMI observations 
of tilted systems indicate that dislocation between Tb ]9 and Tbs 5 can be as far as 100 km if there 
is a strong wind shear. Such dislocation not only poses a problem to rainfall retrieval algorithms 



that use only scattering information, but also causes large uncertainty in rainfall retrieval from 
multichannel retrieval algorithms. This suggests that combined radar and radiometer data is 
needed to reduce the effect of tilt and to improve surface rainfall retrieval. 



1. Introduction 


Measurement of tropical rainfall is very important to understand the hydrological cycle and its 
role in the global climate system. Satellite observation of tropical rainfall is necessary due to the 
huge extent of oceans in the tropics. Comparing to visible/infrared observation, passive 
microwave observation provides better physical information for estimating surface rainfall 
because microwaves can penetrate clouds and interacts directly with hydrometeors at lower 
levels. Upwelling passive microwave brightness temperatures are determined by the surface and 
the vertical distributions of ice and liquid contents in clouds. At low microwave frequencies (< 
37 GHz), the brightness temperature responds to emission from rain and cloud liquid and thus 
the observed brightness temperature increases over the radiometrically cold ocean. In contrast, at 
high frequencies (> 37 GHz), the brightness temperature responds to scattering from cloud ice 
and the observed brightness temperature decreases over both ocean and land backgrounds. 
Rainfall retrieval from satellite passive microwave brightness temperatures has been investigated 
by many researchers (Wilheit et al. 1977, Spencer et al. 1986, 1989, Kummerow et al. 1994, and 
Smith, et al. 1 994). 

Tropical precipitation generally occurs in the form of organized mesoscale convective systems 
(MCSs), which can be characterized by two distinguishable components: convective and 
stratiform (C/S) regions (see Zipser, 1977, Leary and Houze, 1979 and Houze, 1989 for the 
detailed descriptions of MCSs). Convective regions range from a few km to about 30 km in scale 
and have strong updrafts and downdrafts while stratiform regions may extend for hundreds of km 
and have relatively weak vertical air motion and light precipitation. Strong ice scattering has been 
observed by passive microwave radiometers to be associated with deep convection (Wilheit et al. 
1982) due to the strong updrafts in convective regions. However, McGaughey, et al. (1996) has 
found significant ice scattering in stratiform precipitation regions. In such cases, wind shear 
causes the cloud top ice particles to shift horizontally away from heavy surface rainfall regions. 
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Therefore, if the significant vertical wind shear is present, 85 GHz brightness temperatures 
(Tb 85 ) that respond to cloud ice will be displaced from brightness temperature of 10 or 19 GHz 
channels (Tbio or Tbi 9 ) that respond primarily to surface rainfall and liquid hydrometeors in 
lower layers of clouds. Such tilted convective systems have been observed by aircraft radar and 
radiometer as well as ground-based radar during the Tropical Ocean and Global Atmosphere 
Coupled-Ocean Atmosphere Response Experiment (TOGA COARE). Figure 1 displays one 
case which was observed by the Airborne Rain Mapping Radar (ARMAR) and the Advanced 
Microwave Precipitation Radiometer (AMPR). 

Satellite passive microwave radiometers, such as the Special Sensor Microwave/Imager (SSM/I) 
(Hollinger, et al. 1987) and the Tropical Rainfall Measurement Mission (TRMM) Microwave 
Imager (TMI) (Kummerow, et al. 1998), observe the atmosphere at a viewing angle. Thus, the 
brightness temperatures for a tilted system viewed from one direction of tilt may be different 
than those observed from the opposite direction due to different optical depth. The 
displacement between Tbss and Tbj 9 and the viewing angle dependent brightness temperatures 
may result in errors of rainfall retrieval in a tilted system. Thus, investigation of microwave 
brightness temperatures of tilted convective systems is important for improving satellite rainfall 
retrieval. 

In this paper, brightness temperatures of tilted systems at different channels and different 
viewing directions are investigated. Brightness temperatures are computed for a simple tilted 
cloud model and a cloud resolving model. To account for the effect of tilt, the 3D backward 
Monte Carlo transfer developed by Roberti (1994) is used and will be described in the section 2. 
Brightness temperatures of tilted systems calculated from cloud models will be presented and 
discussed in section 3. The successful launch and operation of the TRMM satellite provides us 
with an opportunity to examine tilted convective systems using collocated radar and radiometer 
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data. In section 4, satellite observations of tilted convective systems from Precipitation Radar 
(PR), TMI and SSM/I are presented. Conclusions and summary' are given in section 5. 


2. Radiative transfer model 

Although 1 -dimensional (1-D) plane parallel radiative transfer models are widely used in 
microwave spectral regions (Wilheit, 1977, Kummerow, 1993), 3-dimensional (3-D) radiative 
transfer models have been developed for horizontally and vertically finite clouds. Weinman and 
Davies (1978) first developed such a model, but they assumed vertical homogeneity. Kummerow 
and Weinman (1988) continued this study to include both liquid and ice layers and found that 
footprint-averaged brightness temperatures from finite clouds deviated considerably from the 
plane-parallel approximation. Haferman et al. (1993) used the discrete ordinates method to 
develop a 3-D radiative transfer model and reported that the nadir brightness temperatures from a 
finite cloud over land were lower than the nadir brightness temperature from a horizontally 
infinite cloud due to net leakage of radiation through the sidewalls of the cloud. Using a reverse 
Monte Carlo model. Petty (1994) described the reflected image of a finite cloud by the surface. 
Roberti et al. (1994) developed a 3-D backward Monte Carlo radiative transfer model and 
compared 1-D versus 3-D microwave radiative transfer. They observed that the brightness 
temperature pattern from 3-D microwave radiative transfer is shifted and enlarged when the 
viewing angle is nonzero. 

The so-called 3-D effects, which is difference of brightness temperature over a finite cloud and a 
horizontally infinite cloud, can be attributed to geometric and physical problems (see Haferman, 
et al. 1996). Geometrically, nonzero view angles lead to a larger brightness temperature pattern 
than rainfall pattern and a horizontal shift between brightness temperature field and surface rain 
field. Physically, radiation from the side walls of a finite cloud and the reflected image of cloud 
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by the surface distorts the brightness temperature pattern from the surface rain field. The 3-D 
effects are also reported by Liu (1996) and Bauer et al. (1998). 

To compute microwave brightness temperatures for a tilted system, a 3-D radiative transfer 
model must be used to account for the effect of the tilt. In this paper, we use the 3-D backward 
Monte Carlo radiative transfer model developed by Roberti et al. (1994). The backward Monte 
Carlo method tracks photons that are received by airborne or satellite-borne radiometers back to 
their sources through the medium, following probabilistic interaction laws which are sampled by 
the selection of numbers from a quasi-random sequence. The photons are either absorbed or 
scattered at an interaction point. If a photon is absorbed, it is considered as being emitted with 
the temperature at the point of absorption. When a photon is scattered, a new direction of travel 
is determined by the phase function, which is assumed to be a Henyey-Greenstein phase 
function having a given asymmetry factor. If a photon collides with the surface, scattering or 
absorption occurs depending on the surface emissivity and the random number. If a photon 
escapes from the upper boundary of the cloud, then it is considered as being emitted with the 
cosmic background temperature of 2.7 K. 

Emission and absorption due to atmospheric gases (water vapor and oxygen) and hydrometeors, 
multiple scattering by hydrometeors, surface emission, and cosmic background radiation are all 
accounted for in this model. A more detailed description of the 3-D backward Monte Carlo 
model can be found in Roberti et al. (1994). The extinction and absorption coefficients of 
atmospheric gases (water vapor and oxygen) are computed using the formulas of Ulaby et al. 
(1981). Both liquid and ice are assumed to follow Marshall-Palmer drop size distributions. 
Radiative properties of atmospheric hydrometeors are computed using Mie theory. 


3. Microwave brightness temperatures of tilted cloud models 
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To examine microwave brightness temperatures of tilted systems, the 3-D backward Monte 
Carlo radiative transfer model is applied to two cloud models. One is a simple finite cloud model 
and the other is the Goddard Cumulus Ensemble model (GCE), which is a cloud microphysical 
model developed primarily by Tao and Simpson (1993). Brightness temperatures at 19 and 85 
GHz in horizontal polarization are computed at different viewing angles. For convenience, they 
will be referred as Tbj 9 and Tbss- Results and discussion are given in this section. 

a. Simple tilted cloud model 

A simple finite cloud model was developed by Haferman et al. (1996) and is depicted in Fig.2a, 
where a 5 1 km x 5 1 km x 0.5 km thick cuboidal cloud is centered over a 5 km x 5 km x 4.5 km rain 
shaft and a 5 km x 5 km x 1 km ice layer is centered over the cloud layer. 5 km x 5 km is chosen 
because the TMI has such a resolution for 85 GHz. Following the Wilheit et al. (1977), the cloud 
model assumes that the relative humidity increases linearly from 80% at the surface to 100% at 
the freezing level of 5 km. The temperature of the atmosphere increases from 0 K at the freezing 
level to 305 K at surface with a lapse rate of 6.5 K/km, and the wind speed at the surface is set to 
0 m/s. In order to account for the effect of tilt, we modified this model by assuming that the 
central rain shaft is tilted 45°. Figure 2b and 2c display side views of the vertical and tilted 
systems. 

For the radiative transfer calculations, the cloud model uses 51 x 51 x 12 grid points with grid 
resolutions of 1 km horizontally and 0.5 km vertically respectively. The top boundary of the 
cloud is considered as cosmic background with a microwave brightness temperature equal to 2.7 
K. The lower boundary is the surface which is assumed to be specular. For the tilt of 45°, the 
rain column will shift one grid point horizontally for every two vertical layers. 
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The 3-D backward Monte Carlo radiative transfer model is applied to the tilted cloud systems. 
Brightness temperatures for 19 and 85 GHz at the incident angle of 53° are computed for a rain 
rate of 20 mm/hr. Figure 3 presents the 19 GHz brightness temperature at the azimuth angles of 
0° (solid line, looking from top right to bottom left) and 1 80° (dash line, looking from top left to 
bottom right). It can be seen that the brightness temperature pattern is larger than surface rain 
field. This is due to the tilt that results in an increase of the radiated area viewed by satellite, and 
also because radiation leaked from the sidewall of the cloud. To better understand these 3-D 
effects, the tilted rain column is overlaid on Fig. 3. 

There are two differences for 19 GHz brightness temperatures when the tilted system is viewed 
from right and left ( i/r = 0° and 1 80°). The first difference is that brightness temperature patterns 
have been shifted away from the surface rainfall fields in opposite directions. As discussed in 
Haferman et al. (1996), brightness temperatures of a central rain shaft observed by a satellite for 
an off-nadir angle (0) is shifted horizontally from the surface rain field by a distance about 
Htan 0 , where H is the mean height of the weighting function defined by Wu and Weinman 
(1984). The discussion of this geometric shift is beyond the scope of this paper. The second 
difference is that the Tb for y/ = 180° is less than the Tb for y/ = 0°. This is resulted from the 
different radiative transfer paths which one passes the cloud ice layer ( yf = 1 80°) while the other 
does not ( y/ = 0°). The ice scattering creates a difference between the two viewing directions 
more than 20 K (see A. and A+ in Fig. 3). 

Figure 4 shows the brightness temperatures of the tilted cloud for 85 GHz. The tilted rain 
column is also overlaid on the figure. Unlike the 19 GHz channel that responds to liquid 
hydrometers at lower layers in the cloud, the 85 GHz channel responds primarily to the ice on 
the top of a cloud. Thus, the 85 GHz brightness temperatures for azimuth angles of 0° and 180° 
depress at the same location where there is ice. However, the low Tbss is shifted from the 
surface rainfall field due to tilt. Besides depression, the Tbss is warmer than the background at 
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the satellite viewing sides (see A+ and A- in Fig. 4). This feature, resulting from the “reflected 
image”, has been discovered by Olson (personal notes in 1989?) and reported by Petty (1994). It 
can be explained by the fact that the emission from liquid water at lower layers is reflected by the 
surface. This phenomena exists in all 85 GHz images shown later. For the azimuth angel i// = 

0°, the emission directly from liquid water at lower layers can also be observed by the satellite. 
The saw-like curve of Tbgf at A+ is probably due to the discontinuity of the tilted cloud layers. 
For the azimuth angel t \f = 180°, the increase of Tbss at A. is caused by radiation from the side 
of the cloud layers and then reflected by the surface. The smoother change of Tbss at B+ and B. 
in Fig. 4 is probably due to radiation leakage at the cloud edge. 

Results from the simple tilted cloud model indicate that brightness temperature varies when 
azimuthal view angle changes. For 19 GHz, the difference could be 30K if ice is present when 
the surface rain field is viewed from the direction of tilt ( y = 180°; see Tb at A. in Fig. 3) and 
viewed from the direction away from tilt ( y/ = 0°; see Tb at A+ in Fig. 3). For 85 GHz, the 

depression of brightness temperatures viewed from both directions are very similar because the 
85 responds to the ice on the top of the cloud. However, the higher Tbss from the “reflected 
image" are located at different sides depending on azimuthal viewing angles. 

b. GCE cloud model 

In section 3a, we have tested a simple tilted cloud model. To better understand the effect of tilt, 
a more realistic dynamic cloud model, the GCE, is used in this section. The model assumes that 
liquid and ice are spherical. The distributions of rain, snow, and graupel (or hail) are taken to be 
inverse exponential with respect to the diameter (D) such that 

N(D) = A r 0 exp(-AD) (2) 
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where N(D) is the number of drops of diameter between D and D+dD per unit volume. No is the 
intercept parameter and A is the slope of the distribution given by: 


A = 




(3) 


V P<h J 

The typical intercept parameters used in the GCE model for rain, snow and graupel are 0.08 cm' 
4, 0.04 cm'4, and 0.04 cm"^, respectively. The density of rain, snow and graupel are 1 g cm'3, 
0.1 g cm' J and 0.4 g cm° respectively. The cloud ice is monodisperse with a diameter of 2x1 0'3 
cm and a density of 0.917 g cm’3. 


A tropical squall line from the development time of 60 minutes to 360 minutes has been 
simulated using the GCE model at a time interval of 60 minutes. The model domain is a 128 x 
128 x 28 grid at a horizontal resolution of 3 km, and vertical resolutions of 0.5 km at the lower 20 
layers and 1 km at the higher 8 layers. The tropical squall line was initialized from an 
environment observed on 22 Feb. 1993 during the TOGA COARE. Readers are referred to Tao 
et al (1993) for additional details. The integrated precipitating liquid and ice of the GCE model at 
the simulation time of 240 minutes is shown in Fig. 5a. Figure 5b shows the vertical profile of 
the precipitating liquid and ice at the cross line of AA’. It can be seen that the squall line is 
upshear-tilted. 


The Monte Carlo radiative transfer model was applied to examine brightness temperature of the 
dynamical cloud model. Figures 6a and 6b present brightness temperatures for 19 GHz at 
azimuth angles of 0° and 180° respectively. One can see the apparent brightness temperature 
difference of the squall line for 19 GHz when viewed from the front ( y/ = 0°) and from the rear 
(y/ - 1 80°). The strong emission of squall line in the leading edge is the combination of radiation 
directly from heavy rain and radiation from rain reflected by the surface (see Fig. 6a), while the 
emission from the leading edge of squall line has been attenuated by passing back through the 
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whole cloud system (see Fig 6b). Brightness temperatures for 85 GHz at azimuth angles of 0° 
and 180° are displayed in Fig. 6c and 6d respectively. Higher than background brightness 
temperatures can clearly be seen along the cloud boundary at the viewing sides. 

Brightness temperatures along the cross line of AA’ in Fig. 5a are displayed in Figs. 7 and 8. 
Figure 7a and 7b compare brightness temperatures at two azimuth angles for 19 GHz and 85 
GHz respectively. Since the shift between the brightness temperature field observed by a 
satellite and the surface rain field due to the 3-D geometric effect may complicate the results 
caused by the tilt of a convective system, we have adjusted the locations of Tb |9 to match the 
surface rain field and Tbss to match the location of ice at the top of cloud. From Fig. 7a, we see 
that different optical paths result in a lower brightness temperatures for y/ — 180°. Due to 
radiation leaked from the side of the cloud, brightness temperatures along the cloud boundary is 
always higher at the viewing side than that viewed from the other side. Figure 7b shows that the 
locations for the strongest scattering of Tbs 5 at y/ = 0° and y/ = 180° differ by about 10 km (see 
B+ and B_). This may be the result of an atmosphere above the ice layer. Emission from the 

“reflected images” can be seen clearly at each viewing side (see A+ and A- in Fig. 7b). 

Figures 8a and 8b compare 19 GHz and 85 GHz brightness temperatures for y/ — 0° and y/ = 
180° respectively. As discussed previously, the 3-D radiative transfer model gives the brightness 
temperature at the top of atmosphere. Therefore, as long as there is a nonzero viewing angle, 
Tbi 9 actually responds to a different vertical profile of hydrometers than the Tbss does. In other 
words, even in a vertically finite cloud system with liquid water in the lower cloud and ice in the 
upper cloud, 3-D computations of high Tbi 9 due to emission and low Tbs 5 due to scattering 
from the same vertical profile will be dislocated. The degree of dislocation depends on the 
thickness of cloud. This makes quantitative examination of displacement between Tbi 9 and Tbg 5 
due to the tilt extremely difficult. But one thing we can conclude is that the tilt of a cloud will 
result in a large displacement between 19 and 85 GHz precipitation signatures. The scattering 
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features are also seen for 19 GHz data, especially for y = 180° (see Fig. 8b). 
c. Different satellite field of view (FOV) 

Sections 3a and 3b give brightness temperatures at the resolutions of cloud models. However, 
what we are interested in is how brightness temperatures change in a satellite field of view 
(FOV). In this section, SSM/I and TMI data are considered. 

The SSM/has seven channels: 19.35 V&H, 22.235 V, 37 V&H, and 85.5 V&H GHz. The TMI 
also has the 19.35, 37 and 85.5 V&H channel pairs, plus two 10.7 GHz channels (V&H), and a 
21.3 GHz water vapor absorption channel instead of the 22.235 GHz channel to avoid saturation 
in the tropics. TMI footprint dimensions range from about 40 km for the 10 GHz channels to 5 
km for the 85.5 GHz channels. The SSM/I has horizontal resolutions from about 48 km for the 
19.35 GHz channels to 13 km for the 85.5 GHz channels. Table 1 lists characteristics of the 
TMI and the SSM/I 

The simulated brightness temperatures in a satellite FOV are computed by convolving the high 
resolution brightness temperature field with the antenna gain function. The antenna gain function 
has been discussed and used by Kummerow et al. (1996) and Olson et al (1996). Figures 9a and 
9b display the 19 GHz brightness temperature from the GCE model at the TMI resolution (~20 
km) for y - 0° and y = 1 80° respectively. Comparing Fig. 9a and 9b with Fig. 6a and 6b, one 
can see that the low resolution TMI data have smoothed features, such as the leading edge of 
squall line, compared to the high resolution images. For the SSM/I resolution, which is presented 
in Fig. 9c and 9d, the squall line is barely apparent when viewed from the rear of the system. 
Although the strong emission of the squall line is still visible at both TMI and SSM/I resolutions 
when viewed from the leading edge, the difference of brightness temperature across the segment 
AA’ for azimuth angles of 0° and 180° is very small (see Fig. 1 la and Fig. 13a). 
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Figures 10a - lOd present Tbss from the GCE model at TMI and SSM/I resolutions for y/ = 0° 
and y/ = 180° respectively. Since the 85 GHz brightness temperature for TMI has a resolution 
of 5km, it almost retains the features at the GCE model’s resolution (see Fig. 10a and 10b). 
Although Tb 85 from SSM/I show the smoothed GCE model features, the higher than background 
Tb 85 due to emission from the side of the cloud are clearly displayed in Fig. 10c and lOd. Figures 
lib and 13b present the Tbss difference across the segment AA’ for azimuth angles of 0° and 
180° at TMI and SSM/I resolutions respectively. Figure 13b indicate that the Tbss difference 
between azimuth angles of 0° and 1 80° for the SSM/I is small. 

Figures 12a and 12b present the displacement between Tbi 9 and Tbss of TMI for azimuth angles 
of 0° and 180° respectively. Figures 14a and 14b are the same as Figs. 12a and 12b but are for 
SSM/I. Note that the dislocation between high Tbi 9 due to emission and low Tbss due to 
scattering is reduced as the resolution degrades. 

The 3-D radiative transfer computation results for Tb (9 and Tbss indicate that the brightness 
temperature of a tilted system varies when viewed from two opposite directions of the tilt. The 
Tb 85 will be displaced from Tbi 9 if wind shear exists. As the resolution degrades to satellite 
FOV, the difference of brightness temperatures viewed from two opposite directions of the tilt 
may be neglected. The dislocation between Tb ]9 and Tbss signatures, however, depends on how 
strong the wind shear is. In other words, the displacement between Tbj 9 and Tbss resulting from 
3-D physical and geometric effects may be neglected for low resolutions of satellite data. But if 
wind shear causes ice particles in the upper cloud to shift more than several satellite footprints 
from the surface rainfall field, the displacement between Tbj 9 and Tbss cannot be ignored (see 
next section). This poses a problem for those rain retrieval algorithms that use only scattering 
information. 
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4. Satellite microwave observations of tilted systems 


Satellite observations of tilted convective systems have been noticed by some researchers. Using 
ground based radar data and infrared data from geostationary satellites, Zipser (1988) and 
Heymsfield and Fulton (1994) found that in propagating squall line systems, the coldest cloud 
tops dislocated from the leading edge to the extensive trailing anvil as the systems evolve. Such 
characteristics were also found during TOGA COARE (Rickenbach et al. 1998). However, 
those observations were limited to locations where there were ground based radars or aircraft field 
experiments. The success of the TRMM project has provided us an opportunity to examine the 
tilted systems globally and frequently, because TRMM carries both the high resolution 
microwave imager, TMI, and the first space-borne Precipitation Radar (PR). The PR is a cross- 
scanning 13.8 GHz radar which has a swath width of about 220 km, about one-third of TMI 
swath width. It can provide the 3-D distribution of rainfall intensity at a horizontal resolution of 
about 4.4 km and vertical resolution of about 250m. The vertical rainfall profile from PR data can 
be used to detect the tilted system. TMI observations of tilted systems at 19 GHz and 85 GHz 
are examined in this section. For comparison, collocated and coincident SSM/I observations are 
also presented. 

The SSM/I data we used can be obtained from the Global Hydrology Resource Center at 
NASA/MSFC (http://ghrc.msfc.nasa.gov). TRMM data are available from NASA/Goddard 
DAAC (http://daac.gsfc.nasa.gov). 

a. Case of January, 05, 1998 

A convective system from 5 January 1998 is depicted in Fig. 15. Figure 15a displays the near 
surface PR reflectivity and Figure 15b shows the cross section of reflectivity along line AA\ 
Figure 15b indicates that the maximum surface rainfall occurs around 30°S and 132.5°W, while a 
large part of ice particles aloft occur over the region west of 1 33°W. Existence of a bright band 
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from 133°W to 133.75°W suggests that the rainfall is strtiform in nature in this region. Thus, ice 
particles do not locate above the convective region, but in the stratiform region. Some were even 
in the precipitation-free region. 

Brightness temperatures for 1 9 and 85 GHz (horizontal polarization) from TMI are presented in 
Figs. 16a and 16b respectively. Note that the warm Tbi 9 in Fig. 16b corresponds to the high 
surface reflectivity shown in Fig. 16a. Fig. 16c shows the Tb] 9 h and Tbssh along the cross line 
AA\ Notice the lowest Tbss along AA’ , which responds to the strongest ice scattering, is about 
1° west away from the highest Tbi 9 . 

Figures 17 a-c are the same as those in Figures 16 a-c but for SSM/I observation of the case. 
Although the brightness temperatures vary more smoothly due to the degradation of spatial 
resolution, the shift between the maximum Tbj 9 and minimum Tbss are still apparent in Fig. 
17c. 

b. Case of August 25, 1998 

Figures 1 8 illustrates the case for the Hurricane Bonnie on 25 August 1998. The near surface PR 
reflectivity and the vertical cross section along line AA’ (from 29°N and 75°W to 32°N and 72°W) 
are shown on Fig. 18a and 18b respectively. A strong convective core appears around 29.75°N 
and 74.25°W, and PR data show apparent attenuation below the height of 3 km. However, a large 
amount of ice particles have moved into the region center over 30.25°N and 73.75°W. 

The brightness temperatures for 19 and 85 GHz channels from both TMI and SSMI are 
presented in Fig. 19 and Fig. 20 respectively. Comparing to Fig. 18, one can see that high 
resolution PR data display a very narrow (about 10 km) rainfall band pattern while TMI data can 
not distinguish such narrow band. Nevertheless, the displacements between 19 GHz and 85 
GHz brightness temperatures along the line AA’ can be seen in both TMI (Fig. 19c) and SSM/I 


13 



data (Fig. 20c). Note that near the longitude 74°W, the cold Tbss occur along the latitude 30°N 
while the warm Tb ]9 appear along 29.75°N. 

Passive microwave observations of the above two cases confirm the results derived from model 
simulations, and show that the dislocation between heavy surface rainfall region and ice particles 
aloft can be as large as 100 km if there is significant wind shear associated with convective 
systems. TMI simulations yield results consistent with observation. Although model 
simulations of SSM/I in section 3 do not show a large displacement between SSM/I Tbi 9 and 
Tb 85 , the satellite observation suggest that even for low resolution SSM/I data, such 
displacement can be appreciable. The displacement of Tbss from surface rainfall is a potential 
problem for single channel 85 GHz surface rainfall retrieval algorithms and also potentially casts 
large retrieval uncertainty to multichannel retrieval algorithms. Such retrieval uncertainty has 
been found in the TRMM operational algorithm 2A12 results for the above two cases (not 
shown). For TRMM operational algorithms, readers are referred to Kummerow et al. (1999). 


5. Summary 

Passive microwave brightness temperatures at 1 9 and 85 GHz of tilted convective systems have 
been examined using data computed from cloud models and measured from satellite radiometers. 
To account for the effect of tilt, a 3-D backward Monte Carlo radiative transfer model has been 
applied to a simple tilted cloud model and to the GCE model to derive brightness temperatures. 
Results from the simple tilted cloud model indicate that brightness temperatures vary when 
azimuthal viewing angle changes. The difference of 19 GHz brightness temperatures over the 
surface rain field is caused by differing optical thicknesses from two opposite viewing directions. 
For 85 GHz, the depression of brightness temperatures viewed from both directions are very 
similar because the 85 GHz channel responds primarily to ice on the top of the cloud. However, 
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the higher than background Tbg 5 due to emission from the side of the cloud are located at 
different sides depending on azimuthal viewing angles. 

Results from the GCE model show similar features. When viewed from the leading edge of a 
squall line, strong emission is apparent due to the radiation leakage from the side of the cloud. 
Viewed from the rear, however, the emission from the heavy surface rainfall field is reflected by 
the ice on the top of cloud so that Tb ]9 is depressed. Due to the existence of atmosphere above 
the ice layer, the 85 GHz depression is shifted about 5 km when viewed from opposite 
directions. Another finding of this study is that warm Tbss, due to emission from the side of 
cloud, exists along the cloud boundary at the viewing side. This phenomena cannot be 
reproduced using a plane-parallel radiative transfer model. 

Because Tbi 9 responds to the low level of clouds and Tbss responds to the upper level of 
clouds, the 3-D radiative transfer model-derived Tbi 9 will be displaced from Tbs 5 if there is a 
non-zero zenith viewing angle. This is true even for a vertical cloud system. But the 
displacement in vertical systems due to the 3-D geometric and physical effects may be ignored 
for low resolution satellite data. Results from section 3c suggest that such displacements are 
reduced dramatically as the resolution degrades. 

The existence of significant wind shear, however, will increase such displacements. Satellite 
observations of tilted systems indicate that dislocation between Tbj 9 and Tbss cannot be 
neglected if there is strong wind shear. Measurements from TMI and SSM/I show that Tbss is 
displaced from Tbi 9 by as much 100 km in the Hurricane Bonnie case. Such dislocation adds to 
uncertainty in rainfall retrieval. 

The investigation of brightness temperatures in tilted systems suggests that care must be taken 
when scattering information from 85 GHz data is used as an indicator of surface rainfall or other 
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index, such as separation of convective and stratiform regions (Anagnostou and Kummerow, 
1996, Hong, et al. 1999). This is particularly true for the higher resolution TMI 85 GHz data. 
Using TMI data only, it may be difficult to detect the direction of tilt and remove the effect of 
tilt in surface rainfall retrieval. Combined with high resolution PR data, however, the effect of tilt 
can be reduced, and the surface rainfall retrieval from satellite data will be improved. 
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Figure Captions 


Figure 1: A tilted convective system observed during TOGA COARE. a) radar reflectivities from 
ARMAR; b) brightness temperatures from AMPR for 10, 19, 37 and 85 GHz. 

Figure 2: A simple tilted cloud model, a) overview of vertical system; b) sideview of the vertical 
rain shaft; and c) sideview of tilted rain shaft. 

Figure 3:19 GHz brightness temperature of the simple tilted cloud model at the azimuth angles 
of 0° (solid line, looking from top right to bottom left) and 180° (dash line, looking from top left 
to bottom right). The tilted rain shaft is overlaid. 

Figure 4: Same as the Fig. 3 except for 85 GHz. 

Figure 5: The GCE cloud model at the simulation time of 240 minutes, a) integrated precipitating 
liquid and ice; b) vertical profile of the precipitating liquid and ice at the cross line of AA’. 

Figure 6: Brightness temperatures from the GCE cloud model, a) 19 GHz, azimuth angle = 0° 
(looking from right side to left); b) 19 GHz, azimuth angle = 180° (looking from left side to right); 
c) 85 GHz, azimuth angle = 0°; and d) 85 GHz, azimuth angle = 180°. 

Figure 7: Comparison of brightness temperatures along the cross line of AA’ at different azimuth 
angles, a) 19 GHz for azimuth angle = 0° (solid line) and azimuth angle = 180° (dash line); b) 85 
GHz for azimuth angle = 0° (solid line) azimuth angle = 180° (dash line). 

Figure 8: Comparison of brightness temperatures along the cross line of AA’ at different 
channels, a) azimuth angle = 0° for 19 GHz (solid line) and 85 GHz (dash line); b) azimuth angle 
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= 180° for 19 GHz (solid line) and 85 GHz (dash line). 


Figure 9: 19 GHz brightness temperatures of the GCE model at satellite radiometer resolution, a) 
TMI, azimuth angle = 0°; b) TMI azimuth angle = 180°; c) SSMI, azimuth angle = 0°; and d) 
SSMI azimuth angle = 1 80° 

Figure 10: Same as the Fig. 9 except for the 85 GHz. 

Figure 1 1 : Same as the Fig. 7 except at TMI resolution 
Figure 12: Same as the Fig. 8 except at TM/I resolution 
Figure 13: Same as the Fig. 7 except at SSMI resolution 
Figure 14: Same as the Fig. 8 except at SSM/I resolution 

Figure 15: Radar reflectivity from PR for the January 5 1998 case, a) near surface reflectivity; 
and b) cross section of reflectivity along line AA\ 

Figure 16 TMI observations of the January 5 1998 case, a) 85 GHz; b) 19 GHz; and c) Tbss 
and Tbj 9 along the cross line of AA\ Lines on a) and b) are PR swath. 

Figure 17 Same as the Fig. 16 except for SSM/I observations. 

Figure 18 Radar reflectivity from PR for the Hurricane Bonnie case on Aug. 25, 1998. a) near 
surface reflectivity; and b) cross section of reflectivity along line AA’. 


21 



Figure 19 TMI observations of the Hurricane Bonnie case on Aug. 25, 1998. a) 85 GHz; b) 19 
GHz; and c) Tbg 5 and Tb ]9 along the cross line of AA\ 

Figure 20 Same as the Fig. 19 except for SSM/I observations. 
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Table 1 . Characteristics of the TMI and the SSM/I 


Instrument 

Frequency (GHz) 

Spatial Resolution (km) 

Polarization (H, V) 


19.35 

69x43 

H, V 


22.235 

50x40 

V 

SSM/I 

37 

37x29 

H, V 


85.5 

15x13 

H, V 


10.7 

63x37 

H, V 


19.35 

30 x 18 

H, V 

TMI 

21.3 

23 x 18 
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37 

16x9 

H, V 


85.5 

7x5 

H, V 













Brightness Temperature (K) Brightness Temperature (K) 





















■S' 






'IP 1 




® r ji 

m 





m- 


■ 


* 

% 




f 















Brightness Temperature (K) Brightness Temperature (K) 


Fig. 8 


a. ¥ = 0 ° (19 & 85 GHz) 



X (km) 


b. 'F = 1 80° (19 & 85 GHz) 






a. TMI (19 GHz, 4* = 0°) k 


BS 



Distance in X (km) 


c. SSM/l(19GHz, ¥ = 0°) K 



Distance in X (km) 


b. TMI (19 GHz, V F= 180°) 


312 - 
259 
206 - 

152 

98 _ 
74 


T T 



i I i J 

125 176 227 278 

Distance in X (km) 


d. SSM/I (19 GHz, ¥ = 180°) 


312 “ 
259 

206 - 

152 - 

98 _ 
74 


T 



125 176 

Distance in X 


227 

(km) 






Distance in Y (km) Distance in Y (km) 


a. TMI (85 GHz, ^ = 0°) 


312 


259 


206 

iljjl ; ' y* 

m 



152 


98 



74 125 176 227 278 

Distance in X (km) 



b. TMI (85 GHz, V= 180°) 

K 

312 


J0U.| 

E 259 

Ji 

250.1 

>- 

i 206 

o 

mm 

200 l 

u 

cc 

mm:. 

m 

Q 152 

m 

150. * 

98 


100. 


74 125 176 227 278 

Distance in X (km) 


c. SSM/I (85 GHz, ¥ = 0° ) k d. SSM/I (85 GHz, ¥ = 180°) k 



Distance in X (km) Distance in X (km) 


Fig. 10 










Brightness Temperature (K) Brightness Temperature (K) 


Fig. 12 







Brightness Temperature (K) Brightness Temperature (K) 


Fig. 13 


a. SSMI 19 GHz pF = 0° & 180°) 



b. SSMI 85 GHz (¥ = 0° & 180°) 




Brightness Temperature (K) Brightness Temperature (K) 


Fig. 14 


a. SSMI 'F = 0° (19 & 85 GHz) 



b. SSMI 'F = 180° (19 & 85 GHz) 
















Tb85 (TMI) 




Tb85 (SSM/I) 


Tb19 (TMI) 


b. Tb19 (SSM/I) 


■ 




i*n 


19 GHz 

85 GHz 

1 1 

i 

Vj 


150 


19 GHz 
85 GHz 

1 

lat: 29.00 

29.75 

30.5 

31.25 

32.00 

lat 29.00 

29.75 31 

Ion: -75.0 

-74.25 

-73.5 

-72.75 

-72.0 

Ion: -75.0 

-74.25 -7t 









